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Abstract 
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. We present a scheme to reconstruct the quantum state of a field prepared inside 

O 1 

^S) ' a lossy cavity at finite temperature. Quantum coherences are normally destroyed by 

o ' 

<^ . the interaction with an environment, but we show that it is possible to recover com- 

O 1 

plete information about the initial state (before interaction with its environment), 
making possible to reconstruct any s-parametrized quasiprobability distribution, in 
Ch ' particular, the Wigner function. 
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1 Introduction 



Recently there have been proposals to reconstruct the quantum state of electromagnetic 
fields inside cavities |T], The reconstruction of non- classical states is a central topic 
in quantum optics and related fields, and there have been a number of proposals to 
achieve it (see for instance ||). In fact, the full reconstruction of nonclassical field states 
H as well as of (motional) states of an ion M have been experimentally accomplished. 
The reconstruction is normally achieved through a finite set of either field homodyne 
measurements, or selective measurements of atomic states [|IJ. 

However, the presence of noise and dissipation has normally destructive effects. In 
fact, the reconstruction schemes themeselves also indicate loss of coherence in quantum 
systems 0. Schemes for compensation of losses have already been proposed || and the 
relation between losses and s-parametrized quasiprobabilities has been pointed out in ref 
0. A method of reconstruction of the Wigner function that takes into account losses (at 
T = 0) has also been presented 0. 

We consider here a single mode high-Q cavity where we suppose that a nonclassical 
field is prepared. The first step of our method consists in driving the generated state 
by a coherent pulse. The reconstruction of the field is done after turning-off the driving 
field, i.e. at a time when the cavity field has interacted with its environment at finite 
temperature. We show that by measuring the density matrix diagonal elements and 
properly weighting them, we can obtain directly the Wigner function even at T ^ 0. We 
should remark that to know a state, one has to have information about all the density 
matrix elements (diagonal and off-diagonal), however, with the method presented here 
(see also 0), it is only necessary to have information about diagonal matrix elements. 

2 Master equation and its solution 

The master equation in the interaction picture for the reduced density operator p relative 
to a driven cavity mode, taking into account cavity losses at non-zero temperature and 
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under the Born-Markov approximation is given by || (in a frame rotating at the field 
frequency u) 



§ = + (i) 

where 

lp = (A + C 2 )p (2) 

with 

Lip = yzapa' — a'ap — pa'aj , L 2 p = — yla [ pa — aa' p — paa' J , (S) 



and 



Kp = - l -[HJ), (4) 



where 

H = ih (a*a — actfj . (5) 

a and are the annihilation and creation operators, 7 the (cavity) decay constant, n is 
the mean number of thermal photons and a the amplitude of the driving field. 

2.1 Displacing the field 

The formal solution to Eq. (0) is given by (see for instance ||) 

p(t) = exp[(K + C)t]p(0). (6) 

It is not difficult to show that Eq. @ can be factorized in the product of two expo- 
nentials, one containing the reservoir (super) operators and the other the interaction fl5|), 
the latter one yielding an effective displacement on the initial field. In order to show this 
we calculate the commutator 



[n,£\p = l%p, (7) 
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which allows the factorization 



p(t) = exp(Ct) exp 



- — (1-e^ 2 ) 
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p(0). 



After driving the initial field during a time t, the resulting field density operator will 
read 



p(t)=e ct pp(0), 



where 



pp(0)=D\P)p(0)D(P), 
with the effective (displacing) amplitude 



(9) 



(10) 



P = -2a- 



1 - e^ 2 



11) 
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2.2 Solution to the Master Equation at finite temperature 

We now obtain the density matrix (|J), by defining |J 

j-p = apa\ J+p = a) pa, J^p = a)ap + pb^a + p, 



(12) 



where the superoperators J_, J + and J3 obey the commutation relations [J_, J + ]p = J^p 
and [J3, J±]p = ±2J±p, which can be written as 

Js 



p(t) = e 2 e 



and we have defined 



n(l-e^) 
n[t) 1 + N t ' 



1 + N t 



3 r s+ i(t)j_ 



MO), 



(13) 



(n + 1)(1 -e" 7 *) 
T+N t 



(14) 



with N t = n{l -e" 7 *). 
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3 The reconstruction method 



We now calculate the diagonal density matrix elements < m\p(t)\m > from ([13]) 



(m\pf){t)\m) = J2 

1 i* JV * fc=0n=0 



k \ m 



n \ n 



[T n (t)r- n [r n+1 (t)] 



k—n 



-n-yt 



[1 + N t , 



2n 



(k\p p (0)\k). 



Multipliying the expression above by powers of the weight function Xs, where 



Xs 



[i+Nt]' 2 + r "^) (f^l ~ r «+l(*) 



(15) 



(16) 



and adding over m we obtain 



FCft .) = E xTHM«)M> = [1 + J v,][r_ x ,r fl(t)] g (*l/W0)l*>, (it) 



m=0 L 1 "'tJL- 1 - A^TlHl |t=0 

where we have used the fact that |L0f 



m=o \ n 

If we now multyply F(/3; s) by the quantity 

2[l + iV t ][l-x^(t)] 

7r(s - 1) 

we finally obtain 

2 ~ /s + l^ fe 



(19) 



W(/?;s) = --7— ttE 73T WM0)l*>, ( 2 °) 



which is the s-parametrized quasiprobability distribution ||TT . 

Therefore, by measuring the diagonal elements of the evolved field density matrix, Eq. 
([TBI), we may obtain complete information on the initial state. 

We should remark that in this case (thermal environment), more information on the 
Pm(t) =< m\p{t)\m > is required than in the zero temperature case. Nevertheless, it is 
possible to find a weight function that allows the reconstruction of the initial field. 



3.1 Measuring the photon distribution 

For completness, we suggest a way to measure the photon number distribution P m (/3, t) of 
Eq. fll7D . It is not difficult to show that the atomic inversion for the case of a three- level 
atom in a cascade configuration with the upper and the lower levels having the same 
parity and satisfaying the two-photon resonance condition is given by (see for instance 

0) 

oo 

W(P;t + r)= £P m (/M)cos([2m + 3]Ar), (21) 

m=0 

where A is the atom-field coupling constant. In order to obtain P m from a family of 
measured population inversions, we invert the Fourier series in Eq. (|2l|) , or 

2A r T max 

P m {f3,t) = — dTW{t + T)cos{\2m + 3}\T). (22) 

7T J 

We need a maximum interaction time T max = n/X much shorter than the cavity decay 
time, which implies we must be in the strong-coupling regime, i.e. A >> 7. 

4 Conclusions 

We have presented a method to reconstruct the Wigner function (and in general any 
quasiprobability distribution) of an initial nonclassical state at times when the field would 
have normally lost its quantum coherence because of the interaction with an environment 
at finite temperature. This is an extension of our previous work where we considered the 
interaction with the environment at zero temperature. The crucial point of our approach is 
the apropriate weighting of the evolved (driven and decayed) photon number distribution. 
Driving the initial field immediately after preparation, is not only useful for covering a 
region in phase space but also makes possible, together with the weight function, Xs, to 
store quantum coherences in the diagonal elements of the time evolved density matrix. 
We have shown here that it is possible to find a weight function that allows to reconstruct 
the initial field even at finite temperature, and this is the main result of our paper. 
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Similar conclusions may be reached by employing the method of generating functions 
13| . The possibility of reconstructing quantum states when the system interacts with an 



environment maybe relevant for applications in quantum computing. Loss of coherence 
for such interactions is likely to occur in such devices, and our method could be used, for 



instance, as a scheme to refresh the state of a quantum computer ]14| in order to minimize 
the destructive action of a hot environment. 
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